Simulation of the two-dimensional Potts model using nonextensive statistics 
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The standard Potts model is investigated in the framework of nonextensive statistical mechanics. 
We performed Monte Carlo simulations on two-dimensional lattices with linear sizes ranging from 16 
to 64 using the Metropolis algorithm, where the classical Boltzmann-Gibbs transition probabilities 
were modified for the nonextensive case. We found that the Potts model undergoes a phase transition 
in the nonextensive scenario. We established the order of the phase transition and we computed the 
critical temperature for different values of the Tsallis entropic index. 
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I. INTRODUCTION 

In the last years nonextensive statistical mechanics was 
successfully applied in the study of various physical sys- 
tems, both from theoretical and experimental point of 
view. Regarding magnetic systems these studies refer to 
the well known Ising model 0, [|J , connected with some 
experimental evidences of nonextensivity in the case of 
CMR manganites 

In the present paper we focus on Monte Carlo sim- 
ulation of the standard Potts model using nonextensive 
statistics. It is known from classical statistical mechanics 
that the Potts model has a richer behavior and it includes 
the Ising model as a particular case. Therefor the study 
of this model in the framework of Tsallis statistics can 
be of interest. In the first part of the paper we review 
some theoretical aspects regarding the Potts model and 
nonextensive statistical mechanics. Section ITVl deals with 
the computational methods used for simulations. In Sec- 
tion[V]we will present and discuss the simulation results. 
Some conclusions and final remarks will be drawn in Sec- 
tion [ 



relation J is a coupling constant and 5{si, Sj) is the Kro- 
necker delta: 



xt \ J 1 if s i = s 3 
if Si Sj 



(3) 



The 2-state Potts model is equivalent to the Ising model, 
the numerical value of the critical temperature in the 
former being one-half of its value in the latter. 

In classical statistical mechanics the critical tempera- 
ture of the r-state Potts model is given by the relation 



ln(l + y/r) 



(4) 



Strictly speaking the two-state Potts model is equivalent 
to the Ising model only if we consider the following ex- 
pression for the hamiltonian: 



W = 53[l-2d(« i ,* i )] 



(5) 
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II. THE POTTS MODEL 

The standard Potts model is a generalization of the 
well known Ising model, however it has a much richer 
phase structure. The two-dimensional r-state Potts 
model is a lattice of spins Sj , which can take the following 
values: 



Si = 1,2, ...,r 
The interaction Hamiltonian is of the form: 

H = — J 22 ^( s ^ s j) 



(1) 



(2) 



where the sum is performed over all pairs of first neigh- 
bours in a square lattice of size L x L. In the above 



In the framework of nonextensive statistical mechan- 
ics the simulation results obtained for the 2-state Potts 
model are identical with those obtained for the Ising 
model (see Ref. 0) only if we choose the Hamiltonian 
given by Eq. ([SJ- This case will be discussed in a forth- 
coming article. 

The ferromagnetic version of the Potts model (J > 
0) undergoes a first order phase transition when r > 4 
and a second order phase transition when r < 4. These 
results are known from classical (i.e. extensive) statistical 
mechanics. In the present paper we will investigate the 
behavior of the standard Potts model described by the 
Hamiltonian ([2]) using nonextensive statistics. 



III. NONEXTENSIVE STATISTICAL 
MECHANICS 
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In 1988 C. Tsallis proposed the following formula 
for entropy, which aims to generalize the classical 
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Boltzmann-Gibbs statistics @, @]: 

1 



Sg — K- 



(6) 



In the above relation pi denotes the probabilities of var- 
ious microstates, k is a positive constant and W repre- 
sents the number of microstates compatible with a given 
macroscopic state of the system. In equation ^ q is a 
real parameter, known as the Tsallis entropic index. For 
q = 1 we recover the classical Boltzmann-Gibbs statis- 
tics. 

Nowadays Tsallis' formula is widely used in the study 
of various physical phenomena. Tsallis statistics became 
of interest in the study of nanoscopic systems after the 
equivalence between Hill's nanothermodynamis and Tsal- 
lis thermostatistics has been proved @. In this context 
we consider that applying Tsallis statistics for spin lat- 
tices and magnetic systems can lead to new developments 
in the study of phase transitions. 

The most important property of the Tsallis entropy is 
its nonextensive nature. If we consider a system which 
consists of two subsystems, A and B, then the entropy 
of the global system is given by the relation [8[ : 



S q (A + B) = S q (A) + S q (B) 



Sq(A)S q (B) (7) 



One can observe that for q = 1 the entropy becomes 
extensive. 

In nonextensive statistics the thermodynamic quantity 
associated with a random variable X is given by the q- 
mean value of the respective variable: 



EiPf 

If we introduce the escort probabilities 



Pi = 



n 



then the q-mean value of X can be written as 



(8) 



(9) 



(10) 



This relation has the same form as the one used in clas- 
sical statistical mechanics, but the probabilities Pi are 
replaced with the escort probabilities defined by means 
of relation ([5]). 

The canonical distribution in Tsallis statistics has the 
following form: 



1 



(1 - q)P{Ei - U q ) 



(11) 



energy constraint, Z q is the generalized canonical parti- 
tion function and U q represents the generalized internal 
energy (defined as the g-mean value of the energy). 

The canonical distribution given by relation pT|) is self- 
referential, so it is not suitable for simulations based on 
the Monte Carlo method. To overcome this problem we 
can use the so-called /3 — ft' transformation [10[ : 



(12) 



Using the above transformation we obtain for the canon- 
ical distribution the following expression: 



Vi = ^ 7 [l-{l-q)P'E i Y- 



(13) 



This form of the canonical distribution is suitable for 
Monte Carlo simulations. 



IV. MONTE CARLO SIMULATION OF THE 
POTTS MODEL USING TSALLIS STATISTICS 

Let us consider a transition from a state of energy Ei 
to another state of energy Ej . If we use the escort prob- 
abilities, the acceptance ratio of the transition is given 
by 



Pi 



-*3 



1 - (1 - q)P'Ej 



1 - (1 - q)0'Ei 



(14) 



Based on the above relation we can generalize the stan- 
dard Metropolis algorithm by replacing the Boltzmann- 
Gibbs acceptance probability with the acceptance ratio 
given by equation (jT4"|) . 

The definition of the physical t emp erature in Tsallis 
statistics is still an open problem [Tl1 - [l3| . Based on ex- 
perimental data in Ref. |3| it was suggested that we can 
consider 1//3' as a temperature scale. This approach was 
used in Ref. |2j by Crokidakis et. al. in the simulation of 
the two-dimensional Ising model. In the present article 
we will follow the same path. 

We performed Monte Carlo simulations on two- 
dimensional lattices with linear sizes of 16,24,32,48 and 
64. For the Tsallis entropic index we considered the fol- 
lowing values: 1, 0.95 and 0.9. As mentioned above, we 
assumed that 1//3' can be used as a temperature scale. 
For every value of ft we performed 10 6 equilibration steps 
and 7.2 • 10 7 production steps. The hardware setup used 
for the simulations consists of 18 Intel Xeon processor 
cores. The software was written in the C programming 
language, based on MPICH2 14] for multiprocessing ca- 
pabilities and SPRNG [15| for random number generation. 



SIMULATION RESULTS 



where pi is the probability corresponding to the energy 
level Ei, (3 is the Lagrange multiplier associated with the 



First of all let us refer to the simulation results ob- 
tained for the 10-state Potts model (r = 10). Figure Q] 
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FIG. 1. Plot of the q-mean value of energy and magnetization as function of temperature in the case of the ten-state Potts 
model, for L — 64 and different values of q. 



shows the plot of the g-mean value of energy and mag- 
netization as function of temperature for L — 64 and 
different values of the Tsallis entropic index. We can ob- 
serve that the phase transition is present even for q =/= 1 
and the critical temperature has a strong dependence on 
q. For the Potts model the magnetization is defined by 
means of the relation: 



Eli M*, i)-i] 

N(r — 1) 



M = 



(15) 



where N represents the number of spins in the lattice. 

In figure [2] we plotted the heat capacity as function of 
temperature for different lattice sizes and q = 0.95. We 
can notice that the peak becomes sharper as the system 
size increases. 
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FIG. 2. Plot of the heat capacity as function of temperature 
in the case of the ten-state Potts model, for q = 0.95 and 
different lattice sizes. 

To evaluate the critical temperature of the transition 
we plotted the fourth order cumulant of the energy as 



function of temperature. 

V E {L) = l- 



(E 4 ) q 

3(£2)2 



(16) 



It is known that Ve{L) presents a minimum at the 
pseudo-critical temperature. Figure [3] shows the plot of 
Ve{L) as function of temperature for q = 0.95 and dif- 
ferent lattice sizes. 
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FIG. 3. Fourth order cumulant of the energy vs. temperature 
in the case of the ten-state Potts model, for q — 0.95 and 
different lattice sizes. 



To establish the order of the phase transition we plot- 
ted the fourth order cumulant of the magnetization (the 
so-called Binder cumulant) as function of temperature 
(see Fig. gj}. 



V M (L) = 1 - 



(M 4 ) q 

m 2 )i 



(17) 



The Binder cumulant exhibits a deep negative minimum 
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FIG. 4. Binder cumulant vs. temperature in the case of the 
ten-state Potts model, for q = 0.95 and L = 64. 



FIG. 6. Binder cumulant vs. temperature in the case of the 
two-state Potts model, for q — 0.95 and L — 64. 



at the pseudo-critical temperature, which indicates a first 
order phase transition [la ]. 

To evaluate the critical temperature in the thermody- 
namic limit, we plotted the pseudo-critical temperature 
as function of 1/L 2 , L being the linear size of the lattice. 
The intercept of the linear fit gives us the critical tem- 
perature in the thermodynamic limit, i.e. for L — >• oo. 
In the case of the ten-state Potts model, for q — 0.95 we 
obtained T c = 0.5909(2). 




in the thermodynamic limit, we plotted the Binder cumu- 
lant as function of temperature for different lattice sizes. 
The curves intersect in one point, which corresponds to 
the critical temperature. 
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FIG. 7. Binder cumulant vs. temperature in the case of the 
two-state Potts model, for q = 0.95 and different lattice sizes. 



FIG. 5. Pseudo-critical temperature vs. 1/L 2 in the case of 
the ten-state Potts model, for q — 0.95. 

For the six-state Potts model the simulation results are 
similar, so we deal with a first order phase transition (for 
q = 1, as well as q =/= 1). 

The two-state Potts model presents a different behav- 
ior. The numerical results obtained by simulations show 
us that this model exhibits a second order phase tran- 
sition. In Fig. [5] we plotted the Binder cumulant as 
function of temperature. The form of the plot indicates 
a continuous (second-order) phase transition. 

To evaluate more accurately the critical temperature 



r = 10 

r = 6 
r = 2 



q = 1.00 

(exact) 

0.701232 

0.807607 

1.134593 



q = 1.00 

0.7012(8) 

0.808(1) 

1.134(2) 



q = 0.95 

0.5909(2) 
0.6864(0) 
0.979(4) 



q = 0.90 

0.4853(2) 

0.569(4) 

0.828(4) 



TABLE I. Critical temperature for different values of r and 

The above table presents the numerical results ob- 
tained for the critical temperature, including all the pre- 
viously discussed cases. One can observe that the critical 
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temperature has a strong dependence on the Tsallis en- 
tropic index q. 

VI. CONCLUSIONS AND FINAL REMARKS 

We performed Monte Carlo simulations of the stan- 
dard Potts model using nonextensive statistics. The mul- 
tiprocessor code used for the simulations was based on 
the Metropolis algorithm, were we replaced the classical 
Boltzmann-Gibbs probabilities with Tsallis acceptance 
ratios. Following the idea presented in Refs. Q and 
[H we used 1//3' as a temperature scale. 

The results show that the two-dimensional Potts model 
undergoes a phase transition in the nonextensive case, 



i.e. for q 1. The order of the phase transition was 
established from the behavior of the Binder cumulant. 
For the ten-state and six-state Potts model we deal with a 
first order phase transition, while for the two-state Potts 
model we have a second order transition. We evaluated 
the critical temperature in the thermodynamic limit for 
different values of the Tsallis entropic index and we found 
that it has a strong dependence on q. 
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